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Abstract 

In this work, we consider a tight binding lattice with two non-Hermitian impurities. The system is described by a 
non-Hermitian generalization of the Aubry Andre model. We show for the first time that there exists topologically 
nontrivial edge states with real spectra in the PT symmetric region. 
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1. Introduction 

The discovery of topological insulators in 2D and 
3D has recently attracted a great deal of attention ]jl}] 
(reference therein). The appearance of the gapless 
edge states within the bulk gap is generally believed 
to be a signature of the topological insulator. Gener¬ 
ally speaking, a topological insulator has gapless ro¬ 
bust edge states while it’s energy spectrum is gapped 
in the bulk. The bulk energy gap closes if a system 
is deformed adiabatically to a topologically nonequiv¬ 
alent system. Remarkably, topological phase is not re¬ 
stricted to two and three dimensional systems. Of spe¬ 
cial importance in the context of localization and topo¬ 
logical phase in ID is the Aubry and Andre (AA) model 
Si IIS S 0 i i 0 El \m 0 • A mapping of the 

ID AA model to topologically-nontrivial 2D quantum 
Hall (QH) system was constructed and used to topo¬ 
logically classify ID quasicrystals SS- Topologi¬ 
cal equivalence between crystal and quasicrystal band 
structures S and between the Fibonacci quasicrystal 
and the Harper model were shown [0. Recently, the 
parameter range for a off-diagonal gapless AA model 
that cannot be mapped onto a QH system was explored 

S. 

The periodic table of topological insulators has been 
constructed only for Hermitian Hamiltonians. This is 
because of the fact that Hermiticity was required in 
quantum mechanics to guarantee the reality of the spec¬ 
trum. However, it was shown in 1998 that Hermitic¬ 
ity requirement is replaced by the analogous condition 
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of PT symmetry, where V and T operators are par¬ 
ity and time reversal operators, respectively ifll . In 
2010, a non-Hermitian system was experimentally re¬ 
alized in an optical system 0. It is well known that a 
PT symmetric non-Hermitian Hamiltonian admits real 
spectrum as long as non-Hermitian degree is below than 
a critical number. Therefore, a natural question arises. 
Does there exist a topologically nontrivial system de¬ 
scribed by non-Hermitian Hamiltonian? This problem 
was considered by some authors [ 16, 1/7,18, 3 20[ 21], 
Hu and Hughes [16] and Esaki et al. [17] studied 
non-Hermitian generalization of topologically insulat¬ 
ing phase at almost the same time. The main find¬ 
ing of these papers is that topological states don’t exist 
in the PT symmetric region because of the existence 
of complex energy eigenvalues. In 0, they consid¬ 
ered Dirac-type Hamiltonians and concluded that the 
appearance of the complex eigenvalues is an indication 
of the non-existence of the topological insulator phase 
in non-Hermitian models. They discussed that finding 
non-Hermitian topological phases without Dirac-type 
Hamiltonians could be possible. Esaki et al. consid¬ 
ered tight binding honeycomb lattice with imaginary 
onsite potentials and non-Hermitian generalizations of 
the Luttinger Hamiltonian and Kane-Mele model 0. 
They found that the zero energy edge states are robust 
against small non-Hermitian perturbation, while these 
states decay because of the imaginary part of eigen¬ 
values. In 2012, Ghosh studied topological phase in 
some non-Hermitian system by changing the metric 
of the Hilbert space i.e. modifying the inner product 
0. Schomerus considered a one dimensional non- 
Hermitian tight binding lattice with staggered tunneling 
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amplitude and showed that the system admits complex 
spectrum 0. Another attempt has recently been made 
to find non-Hermitian Hamiltonian admitting topologi¬ 
cal insulator phase 0. Zhu, Lu and Chen specifically 
considered non-Hermitian Su-Schrieffer-Heeger model 
with two conjugated imaginary potential located at the 
edges of the system |2Q0. They found that zero energy 
modes are unstable and corresponding energy eigenval¬ 
ues are not real in the topologically nontrivial region. 
In all of these papers, no example of topological edge 
states for a non-Hermitian system with a real spectrum 
has been found. The existence of topological phase tran¬ 
sition for PT symmetric non-Hermitian systems is still 
an important question. In the present paper, we con¬ 
sider a non-Hermitian generalization of the off-diagonal 
Aubry Andre model. In contrast to general belief, we 
find topological edge states in the PT symmetric re¬ 
gion, i.e. topological edge states with real spectrum. 
This is the first example in the literature where topolog¬ 
ical states are compatible with PT symmetry. 


2. Model 


We consider ti ght b inding des crip tion of a PT sym¬ 
metric lattice l!22l 23i 24, 25, 2(| 27|]. Let us begin with 
the ID off-diagonal AA model. This is a ID tight bind¬ 
ing model with modulated tunneling parameter. In addi¬ 
tion, suppose two non-Hermitian impurities are inserted 
in the system. The impurities are balanced gain and loss 
materials. We require that the non-Hermitian impuri¬ 
ties are arranged at symmetrical sites with respect to the 
center of the lattice, i.e., particles are injected on the j- 
th site and removed from the (N - j + l)-th site, where 
N is the number of lattice sites. The Hamiltonian reads 


N -1 

+ A cos (2,7i(3n + Q>))a] l a n +1 + h.c. 

n= 1 

+iy(a f j a j - a f N _ j+l a N -j + i) ( 1 ) 

where the parameter y represents non-Hermitian degree, 
t is unmodulated tunneling amplitude, a}, and a„ denote 
the creation and annihilation operators of fermionic par¬ 
ticles on site n, respectively. The first term in the Hamil¬ 
tonian is the kinetic energy from the nearest-neighbor 
tunneling and the last term is the non-Hermitian poten¬ 
tial due to two non-Hermitian impurities. The constant 
A is the strength of the modulation and f3 controls the pe¬ 
riodicity of the modulation. The modulation is periodi¬ 
cal if [A is rational and quasi-periodical if it is irrational. 
As it is usual on the topic of topological insulator in ID, 
the modulation phase <!> is an additional degree of free¬ 
dom. 


Below, we will find energy spectrum of this Hamilto¬ 
nian numerically. We adopt open boundary conditions 
with n — 1 and n = N being the two edge sites. Note 
that the translational invariance of the system is broken 
for open boundary conditions. Before going into the 
details, let us study the PT symmetry of the Hamilto¬ 
nian qualitatively. It is clear that the system has no par¬ 
ity symmetry when the tunneling is modulated quasi- 
periodically. Therefore, the Hamiltonian is not PT 
symmetric when /3 is irrational. Consider now period¬ 
ical modulation of tunneling. For finite number of N, 
the particles on the two ends of the chain may not have 
the same tunneling to their nearest neighbors. Further¬ 
more, particles injected on site j and to be lost from site 
N - j + 1 may feel different tunneling depending on /3 
and N. Therefore, we say that the system with edges 
doesn’t generally possess global PT symmetry in the 
usual sense. For example, for /? = 1/3, the system is in¬ 
variant under usual PT operation when N = 6 but not 
when N = 7 and N — 8. It is well known in the theory of 
PT symmetric quantum theory that a commuting anti- 
linear operator can be found for a non-Hermitian Hamil¬ 
tonian if the spectrum is real |2^ 29j]. The associated 
anti-linear symmetry can be interpreted as PT symme¬ 
try. We will show that our system admits real spectrum 
for some particular cases. 

Our aim is to look for a topologically nontrivial phase 
with real spectrum for our non-Hermitian Hamiltonian. 
Since the PT symmetry is broken for quasi-periodical 
tunneling, we consider periodical tunneling parameter. 
Let us start with (5 — 1/2, where the tunneling alter¬ 
nates at fixed (D between the two values 1 + Tcos(<I>) 
and 1 - A cosLD). Hence, the corresponding spectrum 
has two energy bands. Suppose first that the system is 
Hermitian, y = 0 . In this case, our model is reduced to 
the well known Su-Schreier-Heeger model (SSH) & 
We plot the energy spectrum as a function of the phase 
<I> in the Fig.l. As can be seen, the system has maxi¬ 
mum band gap at O = 0 and the bulk gap closes and 



Figure 1: Energy spectrum for A = 0.4, (3 =1/2 and y = 0 when N = 
49 and N = 50 sites. Topological zero energy edge states appear in 
the whole region of ® when the site number is odd. 
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Figure 2: The real and imaginary parts of energy spectrum for the 
parameters A = 0.4, (A = 1/2, N = 49 sites and y = 0.5. Observe that 
topological zero energy states appear in the unbroken PT symmetric 
phase if j = 2. Note that the spectrum would be real valued when 
7 = 24 if y< 0.423. 


reopens as <J? changes. Remarkably, zero energy states 
appear in the spectrum. It is well known that these zero- 
energy states are localized around the two edges of the 
system 0]. This is a signature of topologically nontriv¬ 
ial phase. The region of nontrivial topological phase 
depends on whether the number of lattice sites is odd 
or even. If N is odd, zero energy edge states appear for 
all <t>. However, this is not case if N is even. The zero 
energy edge state exists at = 0 and varying <I> leads to 
the topological phase transitions that occur at = n/2 
and <t> = 37r/2. The region n/2 < < 3zr/2 corresponds 

to topologically trivial phase. Let us now study the exis¬ 
tence of edge states if the two non-Hermitian impurities 
are introduced into the system. We are looking for topo¬ 
logical edge states with real energy eigenvalues. We 
expect that the effect of impurities on the real part of 
energy spectrum is not appreciable provided N is large 


and y is not much bigger than the tunneling parameter t. 

- To check this idea and study the reality of the spectrum, 
we present numerical calculation for odd and even num- 

m jyer of lattice sites. We see that zero energy states exist 
even in the presence of the non-Hermitian impurities. 
However, we are not interested in zero energy states in 
broken PT symmetric region. We are looking for such 
states in the unbroken PT symmetric region. The key 
point here is the site position of particle injection and 
removal. More precisely, where we inject and remove 
particles has almost nothing to do with the real part of 

“ spectrum but plays a vital role on the reality of the spec¬ 
trum. This can be seen clearly in the Fig.2 (N = 49) 
and Fig.3 (N = 50). Let us start our discussion with the 
case of j = 1, i.e., the particles are injected to the one 
edge and lost from the other edge. In this case, the PT 
symmetric phase is spontaneously broken and the spec- 

_ trum becomes complex in the topologically nontrivial 

m. <i> 

— region whereas the spectrum is real as long as y is be¬ 

low than a critical number in the topologically trivial 
region. As a result, we say that the spectrum is com¬ 
plex in the whole region of the phase 0 < < 2n (in 

the region n/2 < <J) < 3zr/2) when N is an odd (even) 
number. Secondly, let us now suppose that the gain/loss 
occurs at j = / = m/p, where m = 1,2,.. are positive 

^Integers (/ = 2,4,6,.. for (3 - 1/2). It is remarkable 
that the spectrum becomes real in the whole region of 
the phase provided that non-Hermitian degree y is be¬ 
low than a critical number, y c . The critical numbers 
for even and odd values of N are almost the same and 
slightly changes if we increase N. However, it decreases 
with the modulation strength A and vanishes if A > 1. At 
fixed N and A, y, decreases also with increasing m. The 
maximum value of y c occurs when the non-Hermitian 
impurities are on the nearest neighbors of the edges, i.e., 
at j = 2. (y c = 0.56 for A = 0.4 and decreases with in¬ 
creasing /I). As a final step, we numerically confirm that 
the zero energy states are localized around the two edges 
of the lattice. The main finding of this paper is that our 
system has nontrivial topological edge states in the un¬ 
broken PT symmetric region. To this end, we note that 
y c is either very small or zero when j is not equal to /. 
The only exception is the one when the impurities are 
on the two symmetrical neighboring sites with respect 
to the center of the lattice. In this case, y c = 0.43 when 
N = 49 and y c = 0.6 when N = 50. If we change N, the 
critical value y c changes drastically. 

We show that the position of non-Hermitian impurities 
has dramatic effect on the reality of the spectrum for 
P — 1/2. We can extend our previous analysis for other 
values of p - \ /p, where p > 2 is a positive integer. 
In this case, energy spectrum is split into p bands. We 
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stress that topological zero energy edge modes don’t 
exist when p is an odd number. Below we first study 
the reality of the spectrum when p is an odd number 
and then discuss topological phase when p is an even 
number. As an example consider p - 3, where there are 
three bands that are crossed as the phase <b sweeps from 
zero to 2 n. In the presence of impurities, the real part 
of the spectrum is almost left unchanged provided y is 
not much bigger than tunneling amplitude. Compared to 
the case with f5 = 1/2, the impurity lattice sites neces¬ 
sary for the reality of the spectrum are shifted according 
to the formula j = / = m/fi = 3,6,9,..., where the crit¬ 
ical number y c decreases with increasing /. We note 
that the site number N plays also a role on the reality of 
the spectrum. To understand better, suppose that parti¬ 
cles on the one edge of the system feel tunneling with 
amplitude 1 + A cos(<t>). Depending on the number N, 
particles on the other edge of the lattice feel either the 
same tunneling or one of following tunnelings with am¬ 
plitudes 1 + A cos(2;r/3 + d?), 1 + A cos(;r/3 + <T). So, 
there are three different configurations at fixed f and 
we should calculate the spectrum separately for three 
different consecutive numbers. As an example, we take 
N = 48, N — 49 and N - 50 and j = 3. We nu¬ 
merically see that the spectrum is real when N — 50 
unless y > y c = 0.45 (at A = 0.4) and VT symme¬ 
try is broken since the spectrum is complex for N = 48 
and N = 49. Analogously, we expect a similar pic¬ 
ture if N —* N + 3 since the same configuration is ob¬ 
tained. Therefore, the spectrum is real if N equals to 
...,47,50,53,... and y < y c , where the critical num¬ 
ber y c slightly changes with N. Suppose now that p is 
an even number. We are mainly interested in this case 
since topological zero energy states are available. Let us 
study the reality of the spectrum. As an example, con¬ 
sider p = 4, where the impurity lattice sites necessary 
for the reality of the spectrum are given by / = 4,8, ... 
The spectrum becomes real if N — ...,51,55,59,... and 
the critical values y c changes slightly with N at fixed A. 
The critical value y c ~ 0.15 for p = 4 (at A = 0.4 and 
/ = 4 ) and decreases with p. As a result, we say that 
topological zero energy states exist when p is even and 
we show that the corresponding spectrum becomes real 
under certain conditions. 

It is worthwhile to point out that the topological phase 
in the unbroken VT symmetric phase can be broken by 
the next-nearest-neighbortunneling. We introduce next- 
nearest-neighbor tunneling to the off-diagonal VT sym¬ 
metric AA model model. The new Hamiltonian reads 

N—2 

H' = H + t' ^ cilci „+2 + h.c. ( 2 ) 

M=1 


where H is given by (jT} and t' is the next-nearest- 
neighbor tunneling constant. Suppose the next-nearest- 
neighbor tunneling parameter is site independent. With¬ 
out any loss of generality we assume the next-nearest- 
neighbor tunneling is constant. The energy of the edge 
states are not zero anymore and the particle-hole sym¬ 
metry is lost since next-nearest-neighbor tunneling con¬ 
tributes the energy spectrum perturbatively. Introduc¬ 
ing non-Hermitian impurities into the system breaks the 
VT symmetry spontaneously in the presence of next- 
nearest-neighbor tunneling. Although edge states are 
still robust as a direct result of the topological nature, 
complex eigenvalues appear in the system. Fortunately, 
the imaginary part of energy eigenvalues are very small 
in a broad range of y if j = /. For example, for the pa¬ 
rameters t — 1, t' = 0.1, A = 0.4, p = 1/2, j = 2 and 
N = 50, the imaginary part of energy eigenvalues are at 
the order of 1CT 4 when y = 0.2 and 1(T 3 when y = 0.3. 

So far, we have omitted quasi-periodical modulation 
of the tunneling parameter since the system is not VT 
symmetric when /? is irrational. Here, we prove numer¬ 
ically that the spectrum is complex for this case. In the 
absence of non-Hermitian impurities, the energy spec¬ 
trum is broken into fractal set of bands and gaps due 
to the quasi periodical nature of the tunneling param¬ 
eter. As shown in the fig.4, the bands remains almost 
the same but they are connected by the edge states as 
O shifts adiabatically from zero to 2n. The states within 
the band gap are localized states either on the right or the 
left edge of the lattice, while the states within the bands 
are extended. In the presence of non-Hermitian poten¬ 
tial, we numerically check that the real part of the spec¬ 
trum is almost unchanged as long as the non-Hermitian 
degree is weak. However, we see that complex eigenval¬ 
ues appear whenever y is different from zero. We repeat 
our calculations for different values of site number j and 
we see no signature of nontrivial topological phase with 
real spectrum. 

As an another example where quasi-periodicity breaks 
the VT symmetry, let us consider a non-Hermitian gen¬ 
eralization of the diagonal Aubry Andre model 

N -1 N 

H-- Z ^ n ^n+ 1 h.C. + z V cos (2n[3n + <F)a 4 a„ 

n= 1 n— 1 

+iy(a f j a j - a f N _ j+1 a N - j+ 1 ) (3) 

where we suppose /j is irrational. For this Hamiltonian, 
tunneling parameter is constant but the onsite potential 
is quasi-periodic. It is well known that the Hermitian 
model exhibits a localization transition at V = 2 1 when 
the modulation is quasi-periodic. All bulk eigenstates 
are extended for 0 < V < 2t and localized for V > 2 1. 
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The energy spectrum as a function of <1> consists of frac¬ 
tal set of bands and within the band gap edge states 
appear. In the presence of non-Hermitian impurities, 
the spectrum becomes complex. We perform extensive 
numerical calculation but we did not find topologically 
nontrivial phase with real spectrum. 

We have discussed that topological zero energy edge 
states appear in our system when (5 = 1/2,1/4,... 
To understand the topological origin of the edge 
states for /3 = 1/2, let us rewrite the Hamiltonian 
(1) in the Majorana basis, where the annihilation 
operators at site 2 n and 2 n + 1 are defined by 
®2n — @~2ti 1T2/i, d-2n+ I — 7" 2n 4-1 "1" i@~2n+l fzh. The two 

new operators crj and r ; - are two species of Majo¬ 
rana fermions at site j. Therefore the Hamiltonian 
(1) without the non-Hermitian term becomes Hq = 

7,„ I 7~2n7~2/i-hl) J” A_(CT2 H CT2«-1 ~ T"2nT2n I k 

where A+ = -2it( 1 + cos <t>). This Hamiltonian de¬ 
scribes two uncoupled identical Majorana chains and 
the system has Zo_ topological index when A+ > A . The 
Majorana physics is most easily understood in the limit 
A_ = 0. In this case, the system has two unpaired Majo¬ 
rana operators in each uncoupled chains, cn, cr,v and n, 
Tat. So, zero energy Majorana fermions localized at the 
edges appear. However, non-Hermitian interaction term 
includes a coupling term between these unpaired Majo¬ 
rana operators if non-Hermitian impurities are placed at 
the edges, j — 1. Therefore the corresponding eigen¬ 
values become complex valued. If the non-Hermitian 
impurities are placed such that j = 2, then the non- 
Hermitian term of the Hamiltonian (1) in the Majorana 
basis contains a coupling term 1 T 2 T 2 - ctn-iTn-i. Since 
it doesn’t contain unpaired Majorana operators, the sys¬ 
tem admits real energy eigenvalues if the non-Hermitian 
degree is smaller than a critical value. 

To sum up, we have studied the non-Hermitian general¬ 
ization of Aubry Andre model. In this paper, we have 
shown a novel topological aspect for the non-Hermitian 
generalization of the Aubry Andre model. We have 
shown that the system can undergo a topological phase 
transition in the PT symmetric region by varying the 
lattice position of the non-Hermitian impurities. We 
think photonic crystals could be used to observe the 
topological edge states for our model. Our system is 
the first example in the literature and it is worth to look 
for other topologically nontrivial PT symmetric sys¬ 
tems with real spectra. Another important task is to ex¬ 
tend the periodic table of topological insulator to non- 
Hermitian Hamiltonians. 
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Figure 3: The real and imaginary parts of energy spectrum for the pa¬ 
rameters A = 0.4,/? = 1/2, N = 50 sites and y = 0.5. At this value of 
the non-Hermitian degree, topological zero energy edge states appear 
in the unbroken PT symmetric phase only if the non-Hermitian impu¬ 
rities are placed either on the neighbors of the edges, i.e., j = 2 or the 
neighbors of the center of the lattice, i.e. j = 25. If j = 1 or j = 3, 
the spectrum is not real anymore. Contrary to the case with j = 1, 
the spectrum is complex in both topologically trivial and nontrivial 
regions when j = 3. 


N=49 N=50 

E E 



Figure 4: The spectrum for the diagonal AA model for t = 1, V = 0.4 
/? = Vl3 - 3. The site number N changes the particular values of the 
phase where the bands are crossed. 
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